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1. Introduction 

The study of B-meson physics on the lattice is of special phenomenological interest because 
it allows one to obtain constraints on the CKM unitarity triangle. In the standard global unitarity 
triangle fit, the apex of the CKM unitarity triangle is constrained using lattice input for neutral 
B-meson mixing. Experimentally, Bq-Bq mixing is well measured in terms of mass differences 
(oscillation frequencies) Am^; in the standard model it is parameterized by [1] 

Am, = ^^^r]BSomBjlBB^\v:qV,,\\ (1.1) 

where the index q denotes either a <i- or a 5'-quark, m^^ is the mass of the B,-meson, and V,* and Vth 
are CKM matrix elements. The Inami-Lim function, So [2], and the QCD coefficient, Tjg [1], can be 
computed perturbatively, whereas the non-perturbative input is given by fqBs^ : the leptonic decay 
constant fs^ and the B-meson bag parameter Bb^. Computing the ratio Amj/Am^ is particularly 
advantageous because statistical and systematic uncertainties largely cancel and, moreover, the 
ratio of CKM matrix elements becomes accessible 

Am, ^ mg, ^2 Wts\^ (J 2) 

Amrf niBj \Vtd\'^' 

The non-perturbative input is now solely contained in the SU (3) -breaking ratio 




(1.3) 



An alternative way of constraining the CKM triangle has been proposed by Lunghi and Soni 
[3] that uses the CKM matrix elements Vuh or Vd,, thereby avoiding the tension between inclusive 
and exclusive determinations of both Vub (3cJ) and Vch (2c7). This alternative method, however, 
requires precise knowledge of the decay constant /g as well as BR{B — > tv) and Am^. Moreover, 
Lunghi and Soni point out [4, 5] that already the current data may show signs of physics beyond 
the standard model manifesting themselves as deviations of experimental values for sin(2j3) (3.3<7) 
and BR{B — > t/v) (2.8<7) from the Standard Model. They argue that the most likely sources for 
new physics are in Bq mixing and sin(2j3). 

Therefore, it is timely for the lattice community to determine these parameters precisely and 
prove by using statistically and systematically independent setups that all sources of errors are 
well under control. Currently, there are three determinations for t, using 2+1 -flavor gauge field 
configurations: an exploratory study by the RBC-UKQCD collaboration [6] with non-competitive 
errors and two precise determinations by Fermilab-MILC [7, 8] and HPQCD [9], which however 
both rely on the same set of MILC gauge field configurations. The latter two collaborations also use 
the same configurations to obatin the decay constants /g, and /g^ and find good agreement. Here 
we present the first results of our project to determine B-meson parameters using domain-wall light 
quarks and relativistic fe-quarks. 

2. Computational setup 

Our project is based on using the RBC-UKQCD 2+1-flavor domain wall lattices with the 
Iwasaki gauge action for the gluons. Our first results presented in these proceedings are obtained 
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on the coarser a 0.1 Ifm lattices with a spatial volume of 24^ in lattice units [10], whereas the 
full project will also utilize the finer a ^ O.OSfm lattices (32^) [1 1] (for details see Tab. 1). 
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Table 1. Overview of the gauge field ensembles to be used for this project. 

For the heavy quarks we use a relativistic formulation derived from the "Fermilab action" 
[12] in which we tune the three relevant parameters of the action non-perturbatively. As demon- 
straged by Christ, Li and Lin, this requires only one additional experimental input compared to the 
formulation of Fermilab [13]. The relativistic heavy quark (RHQ) action is given by 



5 = I^„ 



;mo + 7oDo ^ + C 



a{Dy 



ICp 



(2.1) 



where the covariant derivative is denoted by D, F^y is the field strength tensor and we need to tune 
the three parameters niQa, cp and l^. Exploratory studies on how to tune these parameters have been 
performed by Li and Peng [14-16]. The idea is to use experimental values for the spin averaged 
meson mass (in = {ms^ + 3mB*)/4) and the hyperfine splitting (A„, = iub* — niB^) as inputs together 
with the constraint from the dispersion relation that the rest mass m\ equal the kinetic mass m2. 
These three quantities are computed for a set of seven trial parameters determined by making an 
initial guess for moa, cp, and and then varying it by a chosen uncertainty iba{„,(,„ 



(2.2) 



(see Fig. 1). We iterate over the parameters {m^a, cp, ^} until we determine the values that repro- 
duce the known experimental measurements for m, A,,,, m\/m2: 
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(2.3) 



with 



Y3-Y2 Y5-Y4 Y,-Ye 



2o, 



and A = J^i — 7 X [moa, cp, Q' . 



(2.4) 
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In (2.4) we use the vectors F, as shorthand notation for [m, A,„, mi/m2]- obtained for the input pa- 
rameters as given in Eq. (2.2), labeled 1-7 from left to right. Eq. (2.3) assumes a linear dependence 
of the meson masses on the parameters of the action; therefore we must be in a linear regime to 
reliably extract {mo<3, cp, (^}. We stop the iteration once the tuned values lie within our variation 
range. 
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Figure 1. Initial guess for the RHQ parameters and their uncertainties. 
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In order to finally compute S-meson decay constants and mixing parameters as precisely as 
possible we repeated the original tuning presented at Lattice 2008 [14] with increased statistics 
and optimized smeared wavefunction parameters used for generating the heavy quark propagators. 
Moreover, performing the tuning on the same configurations we intend to use for computing weak 
matrix elements such as fs will allow for an improved error analysis in which the coiTclations 
among the three RHQ parameters can be fully taken into account. 

Our method for computing the AS = 2 four-quark operators requires the light quarks to be 
generated with a point source and sink, but for the heavy quarks we are free to explore different 
smearing choices. In addition to point sources and sinks, we tried Gaussian smeared sources/sinks 
and also varied the radius of the Gaussian smearing. As a first guess we chose the radius of the 
Gaussian source to be the rms radius of the bb- and cc states [17]. Later we extended the radius 
and found that r^^^ = 0.634fm gives the best signal as can be seen in Fig. 2. 

Using this setup we obtain the RHQ parameters on the three 24^ ensembles as given in Table 2. 
In contrast to earlier work [14] the determination uses only quantities from the heavy-Ught system. 
We expect these values to be close to our final determination. Moreover, we observe that within 
statistical uncertainties there is no dependence on the light sea quark mass (wi^gg). 
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Table 2. Preliminary determination of the RHQ parameters on the three 24 ensembles with a ~ 0. 1 Ifm. 



3. Computation of decay constants 

Using these newly determined RHQ parameters we compute as a first non-trivial test the decay 
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Figure 2. S/-meson effective masses for different fo-quark spatial wavefunctions; in each case the hght quark has a point 
source. 



constants of the unitary (B/) and strange (B^) mesons via the relation 



/mB„ 



(3.1) 



where is the lattice decay amplitude, Z,^ is the renormalization factor. For the meson 
(domain-wall valence quark has mass of physical 5'-quark) we generated data for the set of seven 
different RHQ input parameters. Therefore we can use similar equations to (2.3) and (2.4) in order 
to determine the decay amplitude at the tuned RHQ parameters: 

RHQ 



^RHQ 



niQa 
cp 



(3.2) 



with 



<I>3 — <I>2 — <5>7 — <I>6 

lOmoa ' 2ac/, ' 2af 



and 



= — 7<i) X [moa, Cp, (^]' . (3.3) 



In Fig. 3 we show the results for (^b, obtained on the ensemble with m[^^ = 0.005. To test the 

Unearity with respect to the input parameters {niQa, cp, ^} we used three different sets of seven 

RHQ parameters all centered around the same point. The vertical black line with the gray error 

band indicates the tuned value of niQa, cp or and allows for a simple estimate of the error in 

^B, due to the uncertainty in each of the three parameters. These plots show that the effect of the 

uncertainty in moa and cp is negligible, whereas contributes an error of about 1%. Moreover, we 

see that ^b, depends linearly on moa, cp and in the range of interest. 

Alternatively, one may simply use the tuned RHQ values to compute the needed coiTclation 

functions. We followed this procedure to calculate the B; meson (domain-wall valence quark mass 

equals the light sea quark mass) the decay amplitude 
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Figure 3. Dependence of on the three RHQ parameters {niQa, cp, ^}; results are shown for the Ughtest sea quark 



Atfer renormalizing multiplicatively at 1-loop [18] we obtain the decay constants which 
are given in Tab. 3 and shown in Fig. 4. As mentioned above, the values for fs, are obtained 
by simulating directly at the tuned RHQ parameters, whereas the /g^ values are extracted using 
Eq. (3.2). The statistical errors on fs^ are therefore larger because they take into account the 
statistical uncertainties in the three RHQ parameters. For a better comparison, we increase the 
errors in /g, in the Fig. 4 by the error due to the statistical uncertainty in the RHQ parameters 
estimated in Fig. 3. Moreover, we emphasize that this computation is performed without 0{a) 
improvement. 



4. Conclusion 

We presented our first results computing S-meson decay constants using domain- wall light and 
relativistic Z^-quarks with all parameters of the RHQ action tuned non-perturbatively. Currently, our 
results still need to be ff{a) improved. For fs^ we expect only a mild chiral extrapolation but of 
course need to perform an extrapolation to the continuum and estimate other systematic errors. 
Despite these caveats the small statistical eiTors and the fact that our central values lie in the same 
ballpark as results of other collaborations indicates the promise of our method. 
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Table 3. Preliminary results for the decay constants 
using 1-loop multiplicative renormalization, but with- 
out 0{a) improvement of the axial-current operator. 
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Figure 4. Preliminary decay constants fs, and fs^ 
computed on our three 24^ ensembles. 
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